Quantum Field Theory of Topological Defects 
as Inhomogeneous Condensates 

Massimo Blasone" and Petr Jizba^* 

" Blackett Laboratory, Imperial College, London SW7 2BZ, U.K. 

and Dipartimento di Fisica dell'Universita di Salerno 

and INFN, L 84 100 Salerno, Italy 

CK"' ^ DAMTP, Silver Street, Cambridge, U.K. 

GN . Contents: 

^H . 

Cj ■ 1. Introduction 

O' 

■ I 2. The Haag expansion in Closed-Tinie-Path formalism 

3. Kinks in two-dimensional Xip^ theory 
^ I 4. Vortices in four-dimensional Xip'^ theory 



5. Conclusions 



OO 



:§: 

1. Introduction 

Oh; 

D . Topological defects play an important role in many physical systems ranging from cos- 

'"^ ' mology to condensed matter. Thus they link apparently unrelated areas characterized by 

very different energy and time scales. 
/\ • The issue of spontaneous defect formation during symmetry breaking phase transitions 

j^ ■ has recently attracted much attention [jT|. As originally pointed out by Kibble 0] and more 

recently by Zurek Q], different regions of a system may be unable to correlate during the 
quench time which characterizes the transition and, as a result, some parts of space may 
remain trapped in the original (symmetric) phase, giving rise to topological defects. 

Although the Kibble-Zurek mechanism gives a reasonable estimate of the defect density 
as a function of the quench time (as confirmed, for instance, by recent experiments on 
superfluid Helium [Q), this picture is essentially phenomenological. It is clear that a full 
understanding of the process of defect formation requires a full quantum field theoretical 
formulation of the problem. There has recently been much progress in this direction [^ , and 
here we give some novel results based on the approach which we are currently developing. 
A more systematic account will be presented in two forthcoming papers B|^ . 



*e-mail: M.Blasone@ic.ac.uk, P.Jizba@damtp.cam.ac.uk. To appear in the Proceedings of the 
37th International School of Subnuclear Physics, Erice, 29 Aug. - 7 Sep. 1999. 



We are inspired by the work done in the 70's by Umezawa et ah |^, who showed how 
sohtons can arise in QFT as result of a locahzed (inhoniogeneous) condensation of particles. 
In this picture the extended objects have an inherently quantum origin. The corresponding 
classical soliton solutions are then obtained in the h ^ limit. 

Following these ideas we construct a QFT of topological defects using the Closed-Time- 
Path (CTP) formalism 0. The latter is a vital ingredient of our construction, as it has a 
natural extension to finite temperature and to non-equilibrium situations (which is relevant 
for example in realistic phase transitions). 

The paper is organized as follows: in Section 2, we give a brief account of CTP formalism 
and of some fundamentals of QFT which are essential to what follows. In Section 3, we 
consider a simple model, namely Xip"^ in 2 D (i.e.l + 1), showing how the kink solution can 
be constructed at zero and finite temperature. In Section 4, we discuss the vortex solution 
for AD complex Xip^ theory. We conclude in Section 5. 

2. The Haag expansion in Closed— Time— Path formalism 

Let us consider the dynamics of an Heisenberg field described by 

ijix)=t[H,i;ix)] , II(x)=z[ff,n(x)], (1) 

where 11 is the momentum conjugate to ip and H is the full (renormalized) Hamiltonian in 
the Heisenberg picture. Assuming that the Heisenberg and interaction pictures coincide at 
some time tj, we can write the formal solution of Eqs.(|T]) as ("m" denotes quantities in the 
interaction picture) 

ij{x) = Z'J^A-\t)^,n{x)A{t) , n(x) 
A(t) = e*(*-*0^?n e-*(*-*')^ , U{t2; h) = 

where ?if„ is the interacting Hamiltonian, T is the time ordering and Z^p, Zu are the wave- 
function renormalizations (usually H oc ■?/', and so Z^ = Zji). Eqs.(0) must be understood 
in a weak sense: if not, we would get the canonical commutator between if) and H equal 
to 2Z^(5^(x — y). This would imply Z^ = 1, whilst the Kallen-Lehmann representation 
requires Z^ < 1. The solution of this problem is well known, the Hilbert spaces for ijj 
and ipin are unitarily inequivalent, and the wave function renormalizations Z^ and/or Zn 
then "indicate" how much the unitarity is violated [|T0|. Thus, the choice of the ipin and of 
the associated Hilbert space is not unique: selecting a particular set of ipin corresponds to 
defining a particular physical situation, i.e. initial-time data for the operator equation (p. 
We will use this feature to construct a QFT which "contains" topological defects. 
If derivatives of fields are not present in Cj (and thus TCi = —Cj) we can write 



Zl/'A-\t) H,„(x) A(t) 




(2) 


t2)A-\ti) = TJexp 


Jtl 


} (3) 
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^{x) = Zj U{ti]t)'4>in{x)U ^{ti]t) = Zj Tc <^i„(x)exp 



d X i^j„ 
c 



(4) 



Here C denotes a closed-time (Schwinger) contour (see Fig.l), running from t, to a later 
time tf and back again. Tc denotes the corresponding time-path ordering symbol. 

In the limit tj —>■ — oo, tpm becomes the usual in-(or asymptotic) field. Since tf is 
arbitrary, it is useful to set tj = +00. Eq.(|[) is the so called Haag expansion, for the 
Heisenberg field ip M. 
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FIG. 1. The Closed Time Path 



Generalization of Eq. (HI) to an arbitrary polynomial P oi ip reads 



T{P[^]} = Z^Tc\P[i'^n]exp 



c 



(5) 



3. Kinks in two-dimensional Xip"^ theory 

In this Section we apply the formal considerations developed above to a specific model, 
namely the 2D Xip'^ theory both at zero and finite temperature. Let us consider the following 
Lagrangian (we adopt the Minkowski signature (H — )): 



c = li^,^r-l^^'r-^r. 



(6) 



The Heisenberg equation of the motion for the field ip is 

(d^ + f^^)^{x) = -\tP%x) . (7) 

For /i^ < 0, this model admits at a classical level kink solutions. Defining 

^lJ{x) = V + p{x) with v = {0\^{x)\0) , -2//^ = m^ -/x^ = A^;^ g = V2X , (8) 
we obtain 



1 ^2 1 

r ^/o\2 -"-22 y 4, ^ 3, 

>C = -(a^p) - -m p - —p - -mgp + 



1 

2' 



m 



-mgp [x) -^9 P [x) ■ 



(9) 
(10) 



2 "^ 8 ^ 2 "^ 16A 
(9^ + rn?)p{x) = 

Note that m"^ > 0. The asymptotic field (at t — * — oo) now satisfies 

{d^ + m^)pUx)=0. (11) 

The Haag expansion for the field ip reads (we suppress the renormalization factor Z^) 



ip{x) = V + Tc < Pin{x) exp 



d^y^Liy) 



c 



(12) 



where £[„ = xPfn + '^QP'tn- We now consider the following (canonical) transformation: 



Pin{x) -^ pI^ = Pin{x) + f{x) , (9^ + m^)f{x) = , 



(13) 



with / being a c-number function. It is known |TT[ that such a canonical transformation 
induces an inhomogeneous condensation of pin quanta. 

As the shifted field pf^ fulfils the same asymptotic equations as the unshifted field pin we 
may equally well use it as initial-time data for ip. In this case the Haag expansion reads 



^/(x) = V + Tc\p{nix)exp 



Yj/y^lniy) 



(14) 



The superscript / in il}f{x) and £j„ indicates that we work with asymptotic field pf„. Using 
the (operatorial) Wick's theorem [|l^], we rewrite ([T^) in the following way 



il)' [x) = V + 



i6J{x) 
X : exp 



+ fix) 



expi TT / d^y^t^ 



i I d yJiy)piniy) 
c 



c 



exp 



i5J 
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d'yd'zJiy)Aciy;z)Jiz) 



c 



(15) 



j=o 



where Ac{x;y) = {Tc[pin{x)pin{y)]) ((•••) = (0| . . . |0) where |0) is the vacuum for the 
Pin field). We now take the vacuum expectation value (vev) of Eq.(^), use the relation 
^ [ih] ^i-^l = ^ [^] ^[^]e*^''li^=o' perform the change of variables K ^ K + f and set 
to zero the term with J (there are no derivatives with respect to it). The result is 



{ip^ {x)) = V + exp[/ia] K{x) exp 



h 



(16) 



K=f 



where we have reintroduced h and defined 

''^-\j/"^'y^^^'^y\K(z)5K(y) ' 



; / d'z 
c 
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^—K\z) + mgK\z) 



(17) 



After some manipulations we arrive at the following form for the order parameter 

(^/(x)) = V + C[K]{x)\^^f + D[K]{x)\j,^^ , (18) 



C[K]{x) 



(fyAc{x,y) exp[^a] exp 



c 



D[K]{x) = K{x) exp[^a] exp 



5K{y) 



V 



)f 



(19) 
(20) 



We remark that the order parameter Eq.(|T8D still contains all quantum corrections. 

Classical kinks at zero temperature 

Let us now deal with the Born, or classical, approximation of ([T^)- For this purpose we 
keep only finite parts in C and D terms in the fl — > limit. We then get: 



oo X 



oo _. 

D{h ^ 0)\finite = Resn^o [Dm = K{x) J^ -—^ a^b'' . 



nrn+l 



a'^b 



n=0 



(21) 



(22) 



n=0 



where iGr^x, y) = 6{xq — y^) A(x, y) is the retarded Green's function of the free theory and 
A{x;y) = {[pin{x) , pin{y)]) is the Pauh- Jordan function: 



A(xo,xi;0,0) 



Finally we may write 



(2vr) 



S{k^ - m^) eiko) e 



-ikx 



6{xq-\xi\) Joim \ xl-xj). (23) 



{4{x))=v + J2Qn[K]{^) 



ra=l 



K=f 



(24) 



where Qi = f and {i,j, k = 1,2 . . . and n > 1. An empty sum is zero) 



Qnix) 



(f'yGR{x,y) 



mg^Qi{y)Qj{y) + -g^ ^ Qi{y)Qj{y)Qk{y) 



i+j=n 



i+j+k=n 



(25) 



This recurrence relation belongs to the class of the so called functional equations of Cauchy- 
Marley's type p!^JTB| whose fundamental solution cannot be expressed (apart from a very 
narrow class of kernels) in terms of elementary functions. This means that we cannot hope 
to resolve ( |25| ) in terms of general K (or /) . We can however obtain analytical solutions of 
the classical version of Eq.(^ - the analytical kinks. These may be obtained if we realize 
that the convolution of the 2D retarded Green's function Gr{x) with an exponential is 
proportional to the very same exponential. Thus, if Qi{x) = f{x) is an exponential, Fourier 
non-transformable, solution of Eq.(|l3D, we have Qn{x) oc {Qi{x))"' and 

(26) 

where 7 = (1 — u^)^2 . Plugging this form into ([25| ) we arrive at the following equation 



Qn{x) = AJ^ix) = A„e±'""^(^^-^«") 



A. 






V i+j=n i+j+k= 



AiAjAjs 



(27) 



A solution of (pTj) is An = 2v (^) , with s being a real constant. The order parameter is 

'sf{x)Y 



00 y 

{4{x)) =v + 2vY,[- 

n=l ^ 



2v 



(28) 



Thus, provided f{x) is an exponential solution of the linear equation (p!3|), the solution ( p8| 
fulfils the (classical) Euler-Lagrange equation of motion: 



,2 v./,/ 



{d'+p'Wo{x)) = -xm^)r 



(29) 



which is nothing but the vev of Eq.(|^) in the Born approximation. For instance, if we choose 
f(x) = e-"^T(^i-^o«) with s 



ii^Ux)) 



-2ve"^"'"', we obtain the standard kink solution |T4--|T6 
= f th 



Ly7((a;i 



XqU 



(30) 



describing a constantly moving kink of a permanent profile with a center localized at a + uxQ. 

Classical kinks at finite temperature 

Let us now concentrate on the finite temperature case. It is important to understand 
what happens with the topological defects if the system is immersed in a heat bath at 
temperature T. 

The crucial observation at finite temperature is that the operatorial Wick's theorem still 
holds [|T^ and thus Eq.(|T5|) retains its validity provided (<^ . . . ^ denotes thermal average) 



together with 



{0\Tc{pin{x)pin{y))\0) -^ Ac{x,y;T) =<^Tc{pin{x)pin{y)) ^ 



(31) 



. : ^ A^(. . .). The thermal normal ordering A^(. . .) is defined is such a 
way [|l^ that <^ A^(. . .) ^= 0. This is of a great importance as all the previous formal 
considerations go through also for finite T. 

At finite temperature the question of h appearance is more delicate than in the zero- 
temperature case. The whole complication is hidden in the thermal propagator Ac(x, y; T). 
To understand this, let us make h explicit. The free thermal propagator in spectral or Mills's 
representation [|l^,|l^ then reads 



i/\c{.x,y]T) = h 



£k 



e '^^^-y'^ p{k)[ec{xo - yo) + h{hk,/T)] 
iAc(a;,y;0) +iA^(a;,y), 



f27r)2 



(32) 



where the spectral density p{k) = (27r)e(/co) 5(/c^ — rn"^) with e^ko) = 6'(fco) — 0{—ko). The 
contour step function 9c{xo — yo) is 1 if yo precedes Xq along the contour C. The Bose- 
Einstein distribution fb{x) = (e^^ — 1)~ . To calculate Res^^Q we need to perform a Laurent 
expansion of /{, around h, i.e. 



l_T_ Ihk^ 
^'^2 hko 12 T ^■ 



(33) 



which converges for h\ko\ < 2nT. Taking a regulator A - ^ 
only soft modes we get (see also ref. pO 

«;^-^(a;) >=t;(r)th 



h ; 



in ko integration, i.e considering 



m{T) 



7( (Xi — a) — XqU 



(34) 



where f^(T) = j (|/i^| — 3A <^: pf„ :^) and 7n{T) = y2\v{T). Thus, at a "critical" temper- 
ature defined by equation: \p^\ = 3A ^: p^„ :^, the kink disappears. 

It is an interesting question to ask whether for higher dimensional systems, this temper- 
ature is related to the critical temperature of the system. 



4. Vortices in four— dimensional Xip^ theory 

In this Section we sketch the treatment of the AD complex Xip'^ theory in presence of 
vortices. For a more thorough discussion see ref. ^. We show how the vortex solution 
can be obtained by "shifting" both the fields of the massless (Goldstone) mode and of the 
massive (unstable) mode: the appearance of a topological charge is controlled by the shift of 
the Goldstone field only. We consider the Ginzburg-Landau type Lagrangian for a charged 
scalar field: 

C = d^^lj^d^'iP - ij,^iIj^iIj - ^iV-Vr ■ (35) 

The (unrenormalized) equations of motion for the Heisenberg operator tp read 

D + ^^ + ^\^(x)A ^(x) = 0. (36) 

We assume symmetry breaking, i.e. /i^ < and {ip) = v = a/— 2/i^/A. We parametrize the 
field as ^jJ{x) = {p{x) + v)e^^^^\ where both the fields p and x ci-re hermitian and have zero 
vev. We also put g = VX and m^ = Xv"^ > 0. The equations of motion for p and x I'ead 

[0-{d^xf + m^]p+^mgp^ + ^gp' = v {d^x? (37) 

d, [{p + vfd'-x] = . (38) 

We now choose as the relevant asymptotic fields the ones described by the equations 

[n - {d^XinY + m^] (p,„ + v) = (39) 

^^[{p^n + Vy^^X^n] = 0, (40) 

i.e. (□ + m^) ipin = , where we have used ipin = {pin + v)e^^^". This means that the 
interaction Lagrangian density is 

C^ = (m^ - h^)^Ij^^Ij - -It/^VI' = '^rn^P - gmp^ - — p' + -mh^ ■ (41) 

The corresponding Haag expansion for the Heisenberg field operator ip is then 



^(x) = (p(x) + v)e'^^''^ = Tc < (Pinix) + t;)e^^'"^"^ exp 



'c 



(42) 



As before, the solutions of the asymptotic equations ( pOf) and (^) are not unique. Indeed, 
we may define the following shifted fields: 

Pin{x) -^ pI,{x) = Pin{x) + f{x) (43) 

Xin{x) -^ xL(a;) = Xin{x) + g{x) , (44) 

satisfying the same equations as the ones for the unshifted fields. The c-number functions 
h = f + V and g solve the coupled equations: 



d,[{{pl) + h'{x))d^g{x)] =0. 
The Haag expansion for ip in terms of the new asymptotic fields reads 



^P^'<^{x) = Tc I (pLix) + v)e'^'«^'^ 



exp 



'c 



(45) 
(46) 

(47) 



By construction the field il)^'^{x) satisfies the same Heisenberg equations as the field il){x). 
A particular choice of solutions of Eqs.(^5D, ( ^6D can lead to the description of a particular 
physical situation, e.g. of a system with topological defects. Observe that in (|47|) the whole 
dependence on g factorises out of Tc- 

Let us now consider the vev of the Heisenberg operator (0), i.e. the order parameter: 

(^^'^(x)) = e'^^^'^Fjix) = e'^(")(0|Tc |(pf„(x) + t;)e'^'"(")exp -t f d'ydj{y) | |0). (48) 



In the classical approximation, the order parameter satisfies 

[n - {d,g{x)Y + m^ + \Ff{x)] Ff{x) = 
d,[F]{x)d^g{x)] = 0. 



(49) 
(50) 

This is nothing but the well known vortex equations |^. A particular solution - a static 
vortex along the third axis - is obtained taking Ff{x) time independent with a radial de- 
pendence only. Eq.(pOD then reduces to the Laplace equation having the polar angle as 
solution: 



9[x) 



n9{x) 



n tan 



X2 

Xi 



(51) 



where the integer n guarantees the single valuedness of the order parameter {%Ij^'^{x)). Ff{r) 
then fulfils the following (static) equation: 






Ff{r) 



XFf{r) . 



(52) 



At this point we can use the fact that the function g appearing in the vortex equations 
(1^), (1^) and in the equations for the shift functions (^5|) , (^61) is the same: we pick up 
the solution of the static vortex equation and plug it into the shift equations (^51), (^Bf ) to 
determine^ the other shift function, /. Let us denote this solution by /. Then the expression 
for the Heisenberg operator in the presence of a static vortex at x = will be 



^In doing this, we use the renormahzed equations, thus the terms with the vev will drop. The 
notation is kept in (|45|), (|4q ) since in general these terms are present (eg. at finite temperature as 
thermal averages). 



^vort.^^^ =Tcl (pin{x) + f{r{x)) + v) e^>^-(^)+*"^(^) exp 



'c 



(53) 



This expression can be used as a starting point for further analysis [0]. Here we only 
consider the simplest approximation, i.e. we completely neglect the radial dependence by 
setting / = 0. Then the following relation holds: 

/(x)~e'"^(")V^(x), (54) 

from which also follows the relation between Green's functions with and without vortices {x 
and y are far from the vortex core): 

G'{x,y) = (0|r(/(x)/t(i/))|0) ~ e*"[''(-)-^(s')] G{x - y) . (55) 

It is interesting to note that a similar relation can be derived for a systems of (charged) 



particles in a Bohm-Aharonov potential. In this case, one has ||22|| : 



Gj^{x,y) =e-^^[^(^)-^(^)]Go(x,i/), (56) 

where the Bohm-Aharonov potential is that induced by a magnetic flux = n times the 
elementary magnetic flux: 



A = VVt , (b ds- A = -2TTn/e . (57) 

In the case of a vortex line of strength n we also have, beside Eq.(^): 

J = nVe , (t ds ■ 3 = 2iTn . (58) 

Thus the "duality" correspondence reads n9{x) <-^ —eQ{x) and J ^^ — eA. 

5. Conclusions 

We have developed a field theoretical (operator) formalism suitable for the description 
of quantum systems containing (topological) defects. The use of the Closed-Time-Path for- 
mulation is crucial in our approach, since it allows to treat systems both at zero temperature 
or in thermal equilibrium as well as systems out of equilibrium. 

We have applied our method to Xip"^ theory in 2D and 4D cases. In 2D case, kink 
solutions were studied at zero and at finite temperature. In 4D case, vortex solution was 
shown to arise from a inhomogeneous condensation of Goldstone modes; an analogy with 
Bohm-Aharonov effect was discussed. 
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